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ABSTRACT 

This article investigates the physics that is responsible for creating the outgoing Poynt- 
ing flux emanating from the ergosphere of a rotating black hole in the limit that the 
magnetic energy density greatly exceeds the plasma rest mass density (magnetically 
dominated limit). The underlying physics is derived from published 3-D simulations 
that obey the general relativistic equations of perfect magnctohydrodynamics (MHD). 
It is found that the majority of the outgoing radial Poynting flux emitted from the 
magnetically dominated regions of the ergosphere is injected into the magnetosphere 
by a source outside of the event horizon. It is concluded that the primary source of 
the Poynting flux is associated with inertial forces in the magnetically dominated re- 
gion proper or in the lateral boundaries of the region. However, the existing numerical 
data does not rule out the possibility that large computational errors are actually the 
primary source of the Poynting flux. 

Key words: black hole physics - methods: numerical. 



There are two known theoretical mechanism for pro- 
' ducing field-aligned outgoing poloidal Poynting flux, S^, 
at the expense of the rotational energy of a black hole 
' in a magnetosphere that is magnetically dominated. There 
' are electrodynamic processes collectively called Blandford- 
Znajek mechanisms in which currents flow virtually paral- 
lel to the proper magnetic field direction (force-free cur- 
rents) throughout the magn etically dominated zon e all th e 
way to the event horizon iBlandford and Znaie^ il977f) : 
iPhinnevI Jl983ft : iThorne et al (llQSeTl . Therefore, these elec- 
trodynamic currents have no source within the magnetically 
dominated black hole magnetosphere. Alternatively, there is 
the GHM (gravitohydromagnetic) dynamo in which large 
relativistic inertia is imparted to the tenuous plasma by 
black hole gravity that in turn creates a region of strong 
cross-field currents (inertial currents), , that provide the 
source of the field-aligned poloidal currents, sup porting 
in an essentially force-fre e outgoing wind (see |Punslj 
i200ll ) or the seminal articles iPunslv and Coronitil lll99CD : 
iPunslvl (llQOlT) ). The force- free electrodynamic current flow 
is defined in terms of the Faraday tensor and the current 
density as F'^'^ Jv = 0. As a consequence, J • E = 0, so 
5*^ in a force-free magnetosphere must be injected from a 
boundary surface. The two types of sources associated with 
these two types of Poynting fluxes are quite distinct in the 
ergopshere: provides the J • E source in Poynting's The- 



orem and the force-free (electrodynamic) component of 5*^ 
emerges from a boundary source at the event horizon. 

Numerical models can be useful tools for understand- 
ing the source of emerging from the ergosphere of a 
black hole magnetosph ere. Some recent 3-D simulations in 
'Dc Villiers ct al'(2003): Hiro se et all j2004) : lDe Villiers etU 
(,2005a b): .Krolik et al (,200^~show emanating from 
magnetically dominated funnels inside of the vortices of 
thick accretion flows. In these simulations, ~ J^(in 
Boyer-Lindquist coordinates which are used throughout the 
following) inside of the ergosphere, since the poloidal fleld 
settle s to a nearly r adial configuration early on in the simu- 
lation [Hiros^^aJ (I2OO4I) . In principle, one can clearly dis- 
tinguish the amount of 5*^ emerging from the ergosphere in 
a simulation that is of electrodynamic origin (as proposed in 
IBlandford and Znaiekl lll977l)) from t he am ount due inertial 
effects (the GHM theorv of lPunslvl j200lM by quantifying 
the relative strengths of emerging from the inner bound- 
ary (the asymptotic space-time near the event horizon) with 
the amount created by sources within the ergosphere. In 
the high spin rate simulations in questions, over 72% of 5*^ 
emerging from the ergopsheric funnel is created outside of 
the inner boundary. 

The simulation that is analy zed, "KDE" is one of a 
family of s o lution s described in |De_Villiers_e^a il i2003D: 
iHirose et aj j2004l) : lDe Villiers et all ^2005alb^ : iKrolik et all 
ll2005r ~The solutions within the family can be distinguished 
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by their spin rate as defined by the ratio of blaclc hole an- 
gular momentum to mass, a/M . The KDE simulation spins 
the most rapidly with, a/M = 0.998. All of the simula- 
tions begin with a torus surrounding the black hole which 
would be stable if not for the ad hoc introduction of loops 
of poloidal magnetic field along the equal pressure contours. 
This destabilizes the torus as the shearing of the loops in 
the differentially rotating plasma creates magnetic torques 
on the plasma that initiates an accretion flow. The angular 
momentum removal by the held is sustained by the inward 
flow of gas that approaches a centrifugal barrier near the 
black hole. This barrier creates an "inner edge" of the accre- 
tion flow that forms a funnel roughly along the gravitational 
equipotential surface. The subsequent accretion of gas is re- 
stricted primarily to the equatorial plane. As the magnetic 
flux loops accrete, the upper portion gets stretched verti- 
cally by gas pressure gradients and electromagnetic forces 
forces away from the hole and the inward part of the loop 
gets severely twisted azimuthally as it approaches the hori- 
zon. The field lines become inextricably tangled. The net 
result is the formation of magnetic field lines that are ap- 
proximately radial, but highly twisted azimuthally, in the 
region near the black hole, restricted to the funnel. 

The KDE solution with a/M = 0.998 has by far the 
most powerful outward Poynting flux in the funnel within 
in the family of simulations. It is the only solution in which 
this energy flux is comparable with the energy flux emerg- 
ing in a mildly relativistic jet created within the funnel wall 
(the boundary or shear layer between accretion disk and the 
funnel). The KDE simulation is the only one that appears 
capable of supplying the high power relativistic jets that are 
observed in kinetically dominated FRII quasars (the kinetic 
luminosity equals or exceeds the bolometric luminosity of 
the accretion flow). As an example in Cygnus A and 3C 82, 
the jet kinetic luminosity is roughly 5-10 times the bolo- 
metric luminosity of the thermal emission f rom the accre- 
tion fl ow (see the supplementary material in lSemenov et a] 
l|20^). Consequently, from an astrophysical perspective the 
KDE simulation is very intriguing. 

Ostensibly, one might think that the spin rate oi a/M = 
0.998 is a benign change from the KDP simulation in 
iDe Villiers et all il2003l) ; iHirose et aJj2004) : lDe Villiers et"ail 
^2005allb^ : lKrolik et all (|200^~^itiio/M = 0.9. A signiflcant 
difference is that the surface area of the equatorial plane 
in the the ergopshere (the active region of the bla ck hole) 
is alm o st 2.5 times larger in the a/M = 0.998 case iPunslvl 
j200ll) : ISemenov et aJ feo04l) . 

The paper is organized as follows. Section 1 introduces 
notation and basic concepts that describe MHD black hole 
magnetospheres. This is not new material, but it is taken 
from existing literature and provides a framework for the 
discussion to follow. The second section presents Foynt- 
ing's theorem in curved space-time in integral form, the 
energy conservation equation for the electromagnetic held. 
The primary result of this section appears in section 2.3: 
over 72% of emerging from the ergopsheric funnel dur- 
ing the course of the simulation is created outside of the 
inner bound ary. This res u lt foll ows directly from the data 
presented in iKrolik et aJ (|200^ and does not require any 
interpretation or manipulation of the data. In order to ver- 
ify that the data sampled is faithful to the simulation in 
iKrolik et all (120051) , we note an independent result from the 



same work in section 2.4: there is also a strong source of 
fleld aligned poloidal angular momentum flux, S'£' in the 
ergosphere. The source of is co-spatial with the source 
of and of similar magnitude, i.e., over 70% of Sf that 
emerges from the ergosphere is created external to the event 
horizon. This means that either the simulation truly has an 
electromagnetic sour ce in the erg o spher e or there are two 
independent plots in IKrolik et aJ ll2005h that have coinci- 
dental errors. In section 3, it is shown that there is a unique 
physical explanation of the source of and S'£ that is 
consistent with Faraday's law averaged over time and az- 
imuth (the sampled data from the simulat ion is averaged 
over time and azimuth in IKrolik et 3 J2005|) l. A physically 
acceptable source must be a cross-fleld poloidal current (in- 
ertial current, since it is associated with strong J x B forces 
in the momentum equation of the plasma). By Ampere's 
law this current supports a toroidal magnetic field that is 
the field component required for a jump in both 5*^ and S'^ 
(averaged over time and azimuth) in the ergosphere. The 
GHM theory of black hole magnetospheres is based on such 
a current. The only thing that would invalidate this straight- 
forward calculation are significant numerical errors in the 
simulation. Section 3 is the only sect i on in this paper that 
manipulates the data in IKrolik et all (|200^ and so it does 
not carry the same weight as the direct results of section 
2. However, the analysis of section 3 is useful (even if it 
turns out to be misleading) for the development of future 
simulations because such an outc ome would indicate gross 
inconsistencies within the data in IKrolik et aJ i2005l) . Sec- 
tion 4 compares the results seen in KDE with other simula- 
tions of rotating black hole magnetospheres. This motivates 
a discussion in the section 5 of the various numerical tech- 
niques and why they tend to either support or not support 
the KDE sim ulation, in particular the d ifferences w ith the 
simula tions in iMcKinnev and Gammie I (|2004); Mc Kinnev I 
i2005h are highlighted. Possible sources of numerical error 
are explored in terms of computational technique and co- 
ordinate systems. There is also a discussion of the need for 
3-D simulations versus 2-D simulations. In the conclusion, 
we stress the need for further simulations to see if this is 
indeed a valid interpretation of the data. Some suggestions 
on how future data should be sampled are offered. Key pa- 
rameters to monitor are mentioned and likely candidates for 
numerical artifacts are noted as a guide to help future efforts 
sort out this complicated turbulent problem. 



1 PHYSICAL QUANTITIES IN 

BOYER-LINDQUIST COORDINATES 

The Kerr metric (that of a rotating uncharged black hole), 
g^v, in Boyer-Lindquist coordinates {r,6,4>,t), is given by a 
line element that is parameterized by the black hole mass, 
M, and the an gular momentum pe r unit mass, a, in ge- 
ometrized units iThorne et "all il98(jl . We use the standard 
definitions, = r'^ + a^ cos^ 6 and A = — 2Mr + a^, 
where A = at the event horizon, = M + V — a^. 
The "active" region of space-time is the ergosphere, where 
black hole energy can be extracted, < r < = 
M + a/ A/2 _ a2 cos^ 6l |Penros3 

The flux of electromagnetic angular momentum along 
the poloidal magnetic fleld direction is the component of the 
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stress-energy tensor, = [l/(47r)]F0a,F°'", in the approxi- 
mation that the field is radial. In steady state, the electro- 
magnetic angular momentum flux per unit poloidal mag- 
netic flux is the toroidal magnet ic fleld d e nsity: —B^ = 
^Z^^"'', where g = -o"^ sin^ IPhinnevI Jl983h : IPuiisivl 
l|200l|) . Similarly, the electromagnetic energy flux along 
the poloidal magnetic field direction, , is the compo- 
nent, Tf"" — \l/(4:n)]FtaF "^ , in the ap proximation that 
the poloidal field is radial Ixhorne et all tlOSfi). In steady 
state, the energy flux per unit poloidal magnetic flux is 
(fl^ /c)B^ , where is th e field line angular velocity 
e3 



iPhinnevI ^1983^ : iPunsI^ J200 j) . Consequently, is useful 
for quantifying the energy and angular momentum fiuxes as 
steady state is approached. From Ampere's law. 
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evaluated at late times (as an approximate steady state is 
reached), one expects simulation data that is averaged over 
azimuth to obey yj—gj" ~ B^^. Therefore, at late times, 
j" is a potential source for the current system that supports 

Even when a system has not reached a time stationary 
state, one can introduce a well-defined notion of that 
becomes the field line angular velocity in the steady state. If 
the field is nearly radial one can simply define the expression, 
Fte = —^IpFe^f,. With this definition, Sl^ is a function of 
space and time and in steady state it become s a constant 
along a perfect MHD fiux tube lPhinnevI Jl983t) . Therefore, 
the EMF across the magnetic field is —ilpFg^ by definition. 



2 THE SOURCE OF ERGOSPHERIC 
POYNTING FLUX IN THE KDE 
SIMULATION 



The simulation "KDE" of iDe Villiers et all J2003D: 
iHirose et all j2004ll : lDe Villiers et all J2005allbh : iKrolik et all 
ll2005fr "is of the most interest since it generates an order 
of magnitud e more than any of the other simulations 
iKrolik et all ll2005h . The magnetically dominated funnel 
spans the latitudes 0° < 8 < 55° at the inner calculational 
boundary, near . Quantifying the magnetic dominance in 
the funnel is the pure Alfven speed, Ua = /(^A-rmfic), 
where n is the proper number density, fi is the specific 
enthalpy of the plasma and B^ is the poloidal field strength: 
within the funnel, 10 < Ul < lO". 



2.1 Poynting's Theorem 

The source of can be determined from the integral 
version of Poynting's Theorem, which is the integral of 
Tt = FtvJ'^ combined with Stokes' theorem. In Boyer- 
Lindquist coordinates, the curved space-time equivalent of 
the " J-E" source of is FtaJ" ■ In steady state, the source 
of is —Q,pFo4,J" in the approximation of a radial field. 
Without the radial approximation, one needs to define a 
poloidal field direction, B^, and a cross-field poloidal EMF, 
E± = —Q,pB^ then the source of is E_lJ^ ■ This can all 
be set up with great mathematically complexity fsee iPunslvl 
i200ll) for sample calculations) and no additional physical 
insight. The reader should remember that the poloidal field 



is not exactly radial in the following and the simple notion 
of J* is used to approximate J^. 

The funnel is highly turbulent and even averaging over 
(j) does not damp the large field strength variations that are 
seen in snapshots of this rotating non-axisymmetric system 
(DeVilliers privative communication 2005). Finding a source 
of Poynting flux in a time snapshot is futile since the field 
strength variability manifests itself as large fluctuations in 
the Poynting flux on the inner calculational boundary near 
the horizon and throughout the ergosphere (Hawley private 
communication 2005). Conversely, over long periods of time 
the total outflow of energy swamps the local random tur- 
bulent behavior se en in any time snapshot. For this reason 
iKrolik et aJ (l2005l) time averaged all the data in order to ex- 
tricate the persistent physics from the strong turbulent fluc- 
tuations that obscure the physics in individual snapshots. In 
order to understand the source, one is interested in where 
the total time integrated Poynting flux that radiates from 
the ergosphere originates. This information can be obtained 
by performing a time integral of the simulation, once the 
large transients have died down. The discrete time average 
and azimuthal average of in the KDE simulation near 
is shown in flgure 1. The time integral of can be ap- 
proximated by multiplying the discrete time average by the 
total time. The contour map in flgure 1 indicates a region of 
strong outgoing in the evacuated funnel, 30° < < 55° , 
r ~ r+. It is clear that suddenly diminishes close to r+ 
at r 1.3A/ - 1.5A/. 

In order to investigate possible source terms for , 
one needs to be explicit about what is plotted in flgure 1. 
The data was averaged over (j> a-nd t (denoted by "<>"). 
Whenever a time lapse of 80 M occurs within the high time 
resolution simulation, data was stored for a time "snapshot" . 
There are 76 of these snapshots that are averaged in flgure 
1, steps 26 through 101. Thus, the discrete time average of 
Poynting's theorem is relevant to figure 1, 



1 r /■ 

i=26 

101 r „ 

i=26 <- •' -I » 

101 r „ 

= [l/(27r^)]^^ / ^aFtcF^ 

+[i/(2^2)]^^ / ^gFt^P'^ 

A — '^a V 2+4 



dV 



' dO d(j) 
dr d(j) 



(2) 



where the four sides of the Gaussian pillbox of integra- 
tion in figure 1 are the four curves labelled "1 -4" and 
dV — ^ —gdrd(j>d6 is the spatial volume element. Periodic 
boundary conditions are imposed at = and cj) = 7r/2. 
Hence, < T^*^ > is equal at both boundaries and does not 
contribute to surface integral on the RHS of (2). 



2.2 Limitations of the Existing Data 

The data that was stored in the run of the KDE simulation 
does not allow us to compute all of the terms in Poynting's 
theorem, (2), in each of the coarse time steps, "i". For ex- 
ample, since at least three consecutive fine resolution time 
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Figure 1. The azimuthally averaged and time averaged (over 75% of the simulation that ends at t = 8080M) Poynting flux from the 
model KDE (a/M = 0.998) of Krolik et al (2005). The figure is a magnification of the inner region of figure 10 of Krolik et al (2005). It 
is an excision of a region, 0° < 9 < 65° , r ~ r-|_ that is a little larger than the ergospheric portion of the magnetically dominated funnel, 
0° < 6 < 55°, r ~ r_|_. The majority of switches-off in a thin layer near r = 1.3M - r = 1.5M (see Krolik et al (2005) for a description 
of the units on the color bar). A Gaussian pillbox, 30° < 9 < 55°, is drawn as a dashed white contour for use in Poynting's Theorem. 
There are 26 grid zones between the inner boundary, r = 1.175M and r = 1.5A/. The plot is provided courtesy of John Hawley. 



steps were not stored at each of the 76 coarse time steps, one 
can not compute the time derivative or the current density 
(from Ampere's law) that appear in Poynting's theorem. 
This highlights the limitation of single snapshot informa- 
tion. A computation of Tt"^, in a single noisy snapshot never 
includes the d/dt{T\) which is likely to dominate in a tur- 
bulent environment, nor can one compute J**. Thus, such a 
calculation is likely to be very misle ading. For this rea son 
the time averaged data presented in iKrolik et all (l2005l) is 
considerably better. Even though J'^ can not be computed, 
one expects that the d/dt{Tt^) term time averages to zero if 
the system is physical (i.e., a finite field energy density in a 
small spatial region should not decay as a persistent source 
of Poynting flux). Even so, we are limited by not being able 
to compute J-E. However, there are conclusions that can be 
reached. First of all figure 1 can be used to compare the total 
radiated from and r^. This is the subject of section 2.3 
and the results do not depend on any analysis of the data, 
but follow directly from the data. Seco ndly, we have time 
averaged information on < > from lKrolik et all ll2005l) 
which allows us to restrict the source terms in (2). Using this 
information in section 3.3, we deduced that there is only one 
source th at is cons istent with the < > and the <Tt^ > 
data from lKrolik e t al (200f2|), < J ■ E >. However, without 
having the actual data, complete with differential informa- 
tion, one can not assess the level to which numerical errors 
invalidate the conclusions of section 3.3. 



2.3 The Horizon is not the Primary Source of 
Poynting Flux 



The total Poynting flux emanating from inner 
ergosphere, r < 1.5M, after the large tran- 
sients have died down, 2000M < t < 8080M, 
is {20/n)Mj2]'^\^J^^FtcF°"-ded<j,. Simi- 

larly, the time integrated 5*^ emanating from in- 
ner boundary (horizon), after the large transients 
have died down, 2000M < t < 8080M, is ^ 
{20/TT)MY,l'^2^J^yr:r^Fto,F''''dedct>. The result of signifi- 
cance from this plot is that y^)°!,,. f,^ y^~gFtaF"^ d9 d(j} > 

3.5 I J2l=26li^/^^t^P'"' '^^d.'^ i' i-s- JS^dt emitted 
from the inner boundary (horizon) during the interval, 
2000 Af < t < 8080M, is much smaller than / dt emitted 
from Ts during the interval, 2000M <t< 8080M. Because 
of the saturation of the dark red color in the plotting 
routine, might be even larger above the switch-off layer 
than indicated in the contour map. Thus, we only have a 
lower bound on the strength of 5*^ above the switch-off 
layer and it is likely that more than 72% of the total 5*^ 
is created within this thin layer. At this point, we can 
conclude that during the interval 2000M < t < 8080M, 
less than 28% of that emerges from the ergospheric 
funnel came from the inner boundary, the rest was created 
external to the boundary. 
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2.4 The Ergospheric Source of Angular 
Momentum Flux 

One could conjecture that without the actual data in hand 
that it might be difficult to assess whether the source seen 
in the false color and the embedded white contours in figure 
1 truly reflects the simulation. Thus, we are compelled to 
find supporting data. First of all, is the external source a 
plotting routine artifact? Secondly, one might ask if there 
is any other supporting information that there is a strong 
electromagnetic source external to the horizon? Fortunately, 
there are a family of line plots of the electromagnetic angu- 
laj momentum flu x as a function of 9 and r in Figure 7 of 
iKrolik et all (|200^. Each line plot represents the averaged 
(in t and 0) electromagnetic angular momentum flux. 



< T/ >= [1/(4^)] [< F^tF'^ > + < F^.gF"^ >] 



(3) 



as a function of 6* on a r=constant surface. Figure 7 shows 
that < > is created predominantly in the ergosphere, 
outside of the horizon. Over 70% of the electromagnetic an- 
gular momentum flux in the funnel is created by sources 
within the ergopshere for 25° < 6 < 55° in a location that 
is cospatial with the source of Poynting flux that is promi- 
nent in flgure 1. Since there are two strong electromagnetic 
sources of similar magnitude in the same location in two in- 
dependent plots of the KDE simulation, it is clear that the 
numerical data r epresents an electro magnetic source. How- 
ever, the data in lKrolik et all i2005ll does not allow one to 
determine if the source is physical or if it is a numerical 
artifact. 



3 POSSIBLE SOURCES OF POYNTING FLUX 

There are three possible source terms for this increase in 5*^ 
when (2) is applied to the pillbox in flgure 1. We discuss 
each of these potential source terms below. 



3.1 Decay of Field Energy 

The second term on the LHS of (2) represents the decay of 
locally stored magnetic fleld energy de/dt, 



dt ^ 



(4) 



There was no data processed from the KDE run that can 
assess the magnitude of this term. In order to get the deriva- 
tive on the RHS of (4) in each of the 76 snapshots requires 
that data be dumped in three consecutive fine resolution 
time steps at each of the 76 coarse resolution time steps 
(clearly 5 consecutive time steps would be superior). How- 
ever, this was not done and the simulation would have to be 
run again with at least three times as much data gathered in 
order to compute this term in (4). Thus, we can not rule out 
the significance of this term. The time averaged field decay 
from a finite spatial region is not a viable physical source of 
long term creation. The existence of a persistent source 
of decaying field strength that does not die off would be an 
indication that the simulation is pathological, spontaneously 
creating local field energy from numerical artifacts. This is 
a serious concern. 



3.2 Latitudinal Poynting Flux 

It is possible that Poynting flux flows in the upper and lower 
sides, 2 and 4, in flgure 1 and gets redirected into a radial 
Poynting flux by some unknown mechanism. This latitudinal 
Poynting flux, 5*^, is defined by the surface terms. 



S^dt; 



101 „ 



^Ftc.F"" dr d(l> 



(5) 



It seems that S would most likely originate in the funnel 
wall, 55° < 8 < 60° in figure 1 as opposed to the evacuated 
polar region of the funnel. Note that there is a very strong 
source of in the funnel wall in figure 1. The region is 
denser than the funnel, Ua ~ 1, hence this putative source 
of would be inertial in character. The conversion of S''^ 
from the funnel wall into in a magnetically dominated 
black hole magnetosphere would be an entirely new source 
of Poynting fiux that has never been postulated before in 
the literature. 



3.3 Inertial Current 

The remaining possible source of in (2) is the inertial 
current term 



J ■ EdVdt 



101 



20 

— M / 

7J- Z — / 

i=26 



Ft4,r + Ftrr + Fter]dv 



(6) 



This is the type of source that exists in GHM |Punsl^j200ll) . 
It is shown below that this is the only source in Poynt- 
ing's theorem that is consistent with the data presented in 
iKrolik et aJ i2005ft if numerical errors are negligible. 

In order to differentiate between the two possible phys- 
ical sources first expand 

r/- = [l/{A^)][Ft^F^^ + FteF"^] . (7) 

Then average the coordinate frame Faraday's law over t and 



< F^o,t > + < Fet.4. > + < Ft4,,B >= 



< F^r.t > 4- < Frt,4, > + < FtS,r >= 



(8) 
(9) 



Notice the simplification that (8) and (9) depend only on 
ordinary d erivatives beca use we are working in a coordi- 
nate frame IPunslvl i200lll . Since < > is approximately 
constant and radial in the the funnel for t > 2000M, 

< F^g^t >~ and < F^r,t >~ in (8) and (9), respectively. 
By the periodic boundary condition on 0, < Fgt,4, >= and 

< Frt,(f, >= in (8) and (9), respectively. Thus, the averaged 
Faraday's law in (8) and (9) reduce to, < Ft^^e >~ and 

< Ftri,,r >~ 0, respectively. These relations can be simply 
integrated to give < F^t >Ri or |< F^g >|> M |< F^t >| . 
Physically, as in fiat space-time rotating axisymmetric mag- 
netospheres, F^t is associated with changes in the poloidal 
magnetic fiux. 

If < F^t > is negligible, as suggested by Faraday's law, 
then the time averaged (7) becomes 



< T/ >~ ^ < FteF""- > , 
47r 



(10) 
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in its application to figure 1. In order to isolate the most 
prominent contributor to the RHS of (10), to the change in 
< > at r ~ 1AM, r ecall the angular momentum source 
in Figure 7 of lKrolik et al (2005). If < F^t > is negligible, 
as suggested by Faraday's law, then the time averaged (3) 
becomes, 



(11) 

has only mild variation 



Since < B-^ >^< Fe4, > 
associated with the geometrical factors arising from < V ■ 
B = >, the major contributor to the source of S£ in the 
ergosphere must be a jump in < F^^ >. Therefore, the only 
physically consistent reconciliation of the increase in both 
< > and < > across the region 1.3M < r < 1.5M 
with (10) and (11) is that the primary source for both is 
a jump in —B^. Ampere's law, (1), averaged over t and 4> 
indicates that J® is the source of jump in — B"^ and therefore 
the ultimate source of both and in the ergosphere. 
In particular, the dominant source term for < > in (2) 
is < Ft^J°' >«< FteJ" >. 



4 COMPARISON WITH OTHER 
SIMULATIONS 

In this section, a comparison is made between the ergo- 
spheric electromagnetic source of KDE and dynamics of the 
ergopshere that are seen in other simulations of rotating 
black holes. 



4.1 3-D Simulations 

The only other 3-D MHD simulation that has been per- 
forme d about a rotating black hole is in iSemenov et all 
l|200^. This solution is flawed in that it assumes a back- 
ground pressure distribution. The simulation follows the ac- 
cretion of an individual flux tube on this background. There 
is no back reaction on the pressure distribution. In this 
limit, it is shown that the relativistic MHD equations for 
individual flux tubes reduce to a nonlinear string equation. 
However, to higher order, the trans-field momentum equa- 
tion is ignored, so the results are not self consistent. The 
simulations are more accurate in the scenario developed in 
ISpruit and Uzdenskvl J2005i") in which patches of nearly ver- 
tical magnetic flux accrete sporadically into to an existing 
central bundle of strong flux. On the positive side, these 
simulations provide much higher resolution than any other 
simulation, as thousands of points are used on each flux tube 
in the ergosphere. This was argued in lSemenov et "SI (l2004l) 
to be required to adequately treat the large gradients in the 
highly twisted and stretched fleld lines. It is therefore the 
only long term simulation in which numerical diffusion does 
not cause the field lines to reconnect in the equatorial plane. 
The artificial reconnection seen in other simulations changes 
the topology of the field and totally changes the large scale 
torques in the system as compressive MHD waves can no 
longer be launched from the potentially powerful compres- 
sive piston in the equatorial plane of the ergopshere. This 
phenomenon produces the most powerful jets from black 
holes (via a GHM process) in purel y analytical treatments 
IPunslvl JioolMPunslv and Coronitil (Il990l) . 



Another positive of the ISemenov et aJ ll2004f) simula- 
tions is that the fiux tube evolution is very clearly displayed. 
As the flux tube accretes inward it is pulled and stretched 
toward the black hole and becomes extremely twisted. In the 
ergosphere, a GHM dynamo (Poynting flux and jet source) 
forms close to the equatorial plane. As the flux tube ac- 
cretes further, the dynamo starts drifting upward and set- 
tles at high latitudes j ust outside the horizon as seen in the 
closeup movies of iSemenov et all JloOdl . The late time loca- 
tion of the dynamo is approximately in the same location as 
the Poynting flux source in figure 1. 



4.2 2-D Simulations 

We review the numerous 2-D simulations starting from the 
most primitive and ending with t he most physically realistic, 
iMcKinnev and Gammie I J2004) : iMcKinnev | ^|2005^ . 



4-2.1 Time Evolution of the Wald Field 

The time evolution of the Wald vacuum field (this is the 
curved space-time equivalent of a uniform magnetic field, 
i.e., the I = l,m = momen t of the vacuum electro- 
magnetic field, see IPunsIvl i200 Jl ) loaded with plasma has 
been of particular interest because it seemed like an as- 
trophysically reasonable poloidal field geome t ry. The first 
simulations of 



Camenzind and Khannal ll2000h : lKoide et all 



<2002h : lKoideH2003ll were short term due to code and com- 
puting limitations. The y produce effects very s imilar to what 
is seen in figure 1 and ISemenov et al i2004l) . A force- free 
electrodynamic (a very tenuous MHD plas ma) simulation 
of the same problem in iKomissarov I i2004t) showed a time 
stationary state that was similar to the above for flux tubes 
that thread the equatorial plane and a Blandford-Znajek 
type solution on the field lines that thread the horizon. The 
field lines that thread the equatorial plane drove an electro- 
dynamic jet by a proc ess very similar to the origin al GHM 
ergospheric disk model IPunslv and Coronitil ilQQ tf). 

Subsequently, a very important MHD si mulation was 
perfor med on the same field configuration in IKomissarov I 
i2005r) that has alt ered our und e rstan ding of this problem. 
The initial state in IKomissarov I (l2005h is a tenuous plasma 
accreting from the outer boundar ies of the simulation grid. 
At firs t J the simulat i on res e mbled ICamenzind and Khannal 
j2000l) : lKorde et al I lliool : iKoide I ll2003l) and looked hke it 
might approach lKomissarov 12003) . However, as mass accu- 
mulated in the equatorial plane, by nearly vertical accretion, 
the density we nt up so as to viola te the force-free approxi- 
mation used in lKomissarov 1 1|2004) . The radial gravitational 
attraction on the dense plasma stretched and dragged the 
once vertical field lines in the equatorial plane radially in- 
ward, toward the horizon. The stretched field lines recon- 
nect in the equatorial plane essentially creating a local split 
monopole near the horizon. The final state is an accretion 
solution on all the field lines that thread the horizon, there 
is no outflow or jet that forms. 

One should note that this result is largely dependent 
on the plasma injection mechanism. In the original model 
of iPunslv and Coronitil il990f) the plasma is created as 
p ositron-electron pairs from a gamma ray cloud as proposed 
in lPhinnevI (Il983l) . As the pairs settle in the equatorial plane 
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they annihilate in a thin disk. The density never gets high 
enough to drag th e field lines inward a s was seen in the ac- 
cretion solut i on of iKomissarov I f2005') and the solution of 
iKomissarov"! (1200411 is more physically representative in this 
scenario. Another interesting initial state is the deposition 
of magnetic flux bundles in the ergosphere that accrete dif- 
ferentially to the plasma as a distinct pha se of a two phase 
accretion flow lSpruit and Uzdenskvl i2005l) . In this scenario, 
the flux tubes are injected into the ergosphere almost ver- 
tically and already possess there own magneto-centrifugal 
outflow. The density in the equatorial plane is therefore low 
and an accretion type solution as in Komissarov (2005) is 
not likely to occur. See the introductory remarks to s ection 
5.3 for a discussion of the simulation in iDe Vilher j (l2005l) 
that shows some of these expected effects. 



i.2.2 2-D Accretion from Tori 

The initia l state in iMcKinnev and Gammiel l)200A : 
iMcKinnevI (120051) is similar to the torus threaded 
with magnetic loop s described in t he int ro duction for 
the s im ulations in iDe Villiers et all i200^|)j_|Hiros^^J| 
J2004I) : iDe Villiers et all j2005allbF r iKrolik et all J2005ir 
There are significant differences, the simulation is 2-D and 
[McKinncv and Gammic (2004); McKinncv (2005) use in- 
going Kerr-Schild coordinates instead of Boyer-Lindquist co- 
ordinates. A different numerical c ode is used as w e ll (see 
iMcKinnev and Gammiel lj2004h : IPe Vilhers et aJ (1200.^ 1 
and references therein). In spite of these differences, the 
late time states are very similar, both evolve to a mag- 
netically dominated funnel threaded by an approximately 
radial field with outgoing Poynting fiux. There is no large 
scale vertical magnetic flux in the equatorial plane of the er- 
gosphere. Both have a powerful inertial outflow along the 
funnel wall. This feature is never de s cribed in m u ch de - 
tafl in IMcKinnev and Gammie I l^2004^ : IMcKinnev I \2m^ . 
so it is hard to compare the funnel wall jet between the 
two sets of simulations, i.e., how much of the outgoing en- 
ergy budget is in the funnel and how much from the funnel 
wall. The main result of interest here is that the simula- 
tions indicate that the force-free, Blandford-Znajek mecha- 
nism drives an out going Poynting flux on field lines in the 

funnel lMcKinnev~ a nd Gammic (200'3). 

The most relevant simulation from IMcKinnev I (|22Q^ 
is when a/M = 0.999. A significant increase in the Poynting 
flux from the funnel, a factor of 5—6, is seen compared to the 
a /M = 0. 9 simulation. Not as dramatic as the increase seen 
in lKrolik" et al ( 2005) . A check of the Poynting flux contours 
has revealed no significant external source of Poynting flux in 
the ergospheric funnel in contrast to the KDE simulation (J. 
McKinney private communication 2005). The implication is 
that the discrepancy arises either from the 2-D versus 3-D 
nature of the simulations, the numerical code choice or the 
the coordinate system chosen. No matter what has caused 
this discrepancy, the major point of interest is whether this 
is a physical difference (i.e., 2-D versus 3-D) or a numerical 
error (coordinates choice or code inadequacies). 



5 POSSIBLE SOURCES OF THE 
DISCREPANCY 

It was noted abov e tha t the simu l ations of 
IMcKinnev and Gammiel lj2004h : IMcKinnev I 1^200^ do 

not support the data presented in figure 1 for the KDE sim- 
ulation. Possible sources of this discrepancy are discussed 
below as a guide to help future numerical efforts that could 
resolve this important physical point. The emphasis is on 
possible sources of numerical errors. 



5.1 Numerical Method 

In section 3.1, the creation of Poynting flux resulting from 
the decay of field energy could not be ruled out as pos- 
sible persistent source. Although this is not physical, it 
could result from numerical artifacts. Th e method used in 
iDe Villiers et al (2003): Hiro se et aJj2004h : lDe Vifliers et al 
l|2005Jb() : lKrolik et al pOOsi) is non-conservative. The field 
is conserved in the sense that V • B = is guaranteed by 
the numerical method and is not a subsidiary condition that 
is imposed. However, global energy is not conserved in each 
time step. This can lead to unphysical jumps in the fluid 
4- velocity which could create unphysical changes in the fleld 
strength through the frozen-in condition, F'^^Uv — 0. This 
means that electromagnetic energy could potentially appear 
and disappear for no physically or mathematically consis- 
tent reason. In regions of fast oscillations and large spatial 
disco ntinuities such a code a l ways generates extra oscilla- 
tions ^^Viliie^^m^^^le^ (l2003h . The extra oscillations 
can allow the field to gain or lose additional amounts of en- 
ergy and energy flux. If this effect is signiflcant there is no 
obvious apparent reason why it should preferentially create 
energy or destroy energy in a time average. Thus, in of itself, 
it would not explain the source in figure 1, unless it is cou- 
pled to another numerical problem (see section 5.2 for exam- 
ple). By contrast the reade r shou l d note that the numerical 
metho d s of .Komissarovl (|2004|) : IMcKinnev and Gammiel 
i2004ll ; IMcKinnev I i2005ri are based on a conservative ap- 
proa ch which is therefore much less prone to this type of 
error iGammie et 11 (l2003h . 

The radial resolution of 26 grid zones between the inner 
boundary and the Poynting flux source in flgure 1 is prob- 
ably adequate. Although, since this feature is suspect more 
grid zones would be preferably in future simulations. More 
of a concern is the 9 resolution. The grid zones were con- 
centrated near the eq uatorial plane since the or iginal intent 
of the simulations in iDe Vifliers et"ail (|2003^: iHirose et 
l)2004h : IpTvifliers et aJ (|2005Jbj) : iKrohk et all (|2005^ was 
to model the accretion disk with high resolution. There 
are 192 zones in the range < 6^ < tt and 26 zones in- 
side the Gaussian pillbox of figure 1. Thus, the resolution 
is not adequate to resolve latitudinal gradients which could 
be very important for understanding the differential version 
of Poynting's theorem. Furthermore, it is important to ad- 
equately resolve the funnel jet wall since this is a highly 
dynamic region that carries most of the outflow energy, yet 
it is just a thin layer. Poor resolution of the funnel wall jet 
could allow energy and mass to diffuse across the fleld lines 
into the funnel in figure 1. 
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5.2 Coordinates 

Another concern is the choice of Boyer-Lindquist coordi- 
nates. The Boyer-Lindquist coordinates have a coordinate 
singularity at the event horizon. This singularity could be 
excised by an inne r boundary as in Do Villicrs ot al ( 2003|): 
iHirose et all j2004 : lDe Villiers et all ll2005Jbl) : iKrolik et all 
ll2005r ~but one still has steep gradients in the metric coef- 
ficients outside the horizon. This has long been recognized 
as a problem for numerical work by gravity wave theorists. 
There is a tendency f or strong reflections to be generated 
at the inner boundary ICampanelli et all (|2003). The sublu- 
minal fast MHD wave speed might be a possible distinction 
from gravity wave physics that allows Boyer-Lindquist co- 
ordinates to be used successfully in MHD. The intention 
was to place the boundary so close to the horizon that all 
information would be flowing inward and reflections would 
not be an issue. For this assumption to hold, the flow must 
be going inward super-magnetosonic (faster than the fast 
MHD wave speed) outside of the inner boundary. How- 
ever, the super-fast nature of the flow was never demon- 
str ated for KDE . A rough order of magnitude calculation 
in iPunslvl (1200111 yields the lapse function, a at the fast 



This puts the fast point 



point, Q = \fKsm6 / ^g^^ « U]^ 
at r Si 1.0651M, at 6 = 30° for example. This is inside of 
the inner calculational boundary at r = 1.175M. The inter- 
action of the flow with the boundary is another potential 
artificial numerical source of Poynting flux. 

In order to resolve this problem, gravity wave 
researchers have gon e to Kerr-Schild coordinates 
ICampanelli et all (1200 ll) . There are no large gradients 
in the metric near the horizon or coordinate singularity. An 
inner computational boundary is placed inside the horizon. 
This method has been employ ed bv .Komissarov . (12003) ; 
iMcKinnev and Gammiel J2004) : iMcKinnev I ^2005^ . The 
only way that information can get from inside or near the 
horizon and react back on the upstream flow is by numerical 
diffusion. 



5.3 2-D versus 3-D 



The c h oice of 3-D bylDe Villiers et all (l2003^ : iHirose et all 
j2004l) : lDe" Villiers et all j2005allbll : lKrolik et aj (|2005^ is cer- 
tainly preferential to 2-D for physical reasons. For, example 
flux tube evolution is vastly different if interchange instabili- 
ties are allowed. These are impossible in 2-D. Unfortunately, 
3-D is much more cumbersome numerically. This is the rea- 
son that the KDE simulation was not sampled adequately 
to determin e the 4-current density. As of t oday, the numer- 
ical code of iDe Villiers and Hawlevl (l2003l) is the only code 
that is efficient enough to produce long term self-consistent 
3-D MHD simulations. One should question how accurate 
the 2-D expedience is at reproducing the physics of 3-D time 
evolution in black hole magnetospheres. The external source 
in figure 1 might be the first indication that 3-D is required 
to capture all of the physical effects. 

The most intriguing question is whether there is a the- 
oretical reason why 2-D can yield different results than 
3-D simulatio ns. For example, the ori ginal ergospheric 
disk model of iPunslv and Coronitil dlQQOt) is axisymmetric 
macroscopically, but microscopically it is 3-D. Buoyant flux 
tubes are created by reconnection at the inner edge of the 



ergospheric disk and recycle back out into the outer ergo- 
sphere by interchange instabilities. The details of this model 
could not exist in 2-D. The 3-D requirement to model the 
physics of vertical flux in the equatorial plane of the er- 
gopshere is sugge sted by new unpublished simulations in 
iDe ViUiersI J2005h . De Vilfiers has placed a Wald field (as 
in section 4.2.1) inside of the torus in the initial state, the 
simulation is o therwise similar to the KDP simulati on with 
a/M = 0.9 in iDe Villiers et all li2003.'l: Hirose et aA 1^2004); 
iDe Villiers et alH2005allbh : lKrolik etaJ l2005l) . A weak Wald 
field gets dragged inward with very little modification to 
what was found for KDP in agreement with ^omissaro^ 
i2005ll . A strong Wald field completely disrupts the torus. 
A moderate strength field can be chosen that is sufficient 
to impede the accretion flow, but not disrupt the disk. As 
accreting flux builds up in the ergosphere, accretion gets 
arrested and large episodic releases of jet energy occur un- 
til the built up magnetic pressure gets relieved (presumably 
by interchange instabilities). The recycling magnetic pres- 
sure resembles the rec ycling of magnetic flux tub es in the 
ergospheric disk model jPunslv and Cororiitil il99(]l) . The cy- 
cle repeats on time scales on the order of the Alfven wave 
crossing time for the ergopshere. This is quite a con t rast t o 
the 2-D simulation of the Wald field in iKomissarov I i2005f) . 
In the 2-D scenario, the vertical flux gets overwhelmed by 
plasma inertia and dragged into the hole. In the 3-D simula- 
tion, the vertical flux impedes the accretion flow and there 
is significant vertical flux that is not pinned inside the hole 
as evidenced by the recycling magnetic pressure in the er- 
gopshere. 

A more general motivation for 3-D is indicated by the 
contrast in the field configurations that give rise to either 
electrodynamic or GHM energy generation in the existing 
simulation literature. 



5.3.1 Electrodynamic Field Configurations 

The ultimate distinction between the electrodynamic solu- 
tions and the GHM solutions is whether the magnetic field 
is organized in the ergosphere so that a tenuous plasma 
can fiow approximately force-free (the electrodynamic con- 
dition). For example, this can occur if a weak field thread- 
ing a dense plasma (i.e., Ua < 1) is deposited near the 
event horizon by an accretion flow. In this scenario, the field 
lines are dragged inward radially and twisted up azimuthally 
by plasma inertia. As plasma is depleted by accretion in- 
ward and magneto-centrifugal and thermal outfiow, the 
flux tube has been essentially pre-configured (twisted and 
stretched) by an initial nonforce-free, inertially dominated 
transient, so that subsequent tenuous accretion can pro- 
ceed approximately force- free. Thi s is the manner i n whic h 
the magnotopshere is built up in |D e_Villiei's et af ("200^); 
iHirosc et al ( 20041: 10 6 Villiers et all J2005a|b() : jKrol i k et lil 
li2(305l) : IMcKinnev and Gammie I (|2004^ : IMcKinnev I BooK 
Alternatively, a strong magnetic field can exist 
near the black hole and can be almost vertical. 



in 'Camenzind and Khanna ' (2000*); "Koide et al " (2002 
jKoide ( 2003); Scmenov ot al (2 004); Komissaro^ (2004 
and the initial state in IKomissarov I (l2005l) . The tenuous 
plasma is restricted from a force-free infiow. However, if an 
inertially dominated transient state (no nforce-free) can form 
as in the dense equatorial condensate of lKomissarov I i2005l) , 
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this will drag the field lines inward radially towards the hori- 
zon and twist them azimuthally in the process. The field has 
then been prepared for future approximately force-free (elec- 
trodynamic) inflow in the ergopshere. 

An initial inertially dominated transient appears to be 
required to modify an accreted field to a force -free configu- 
ration . This is evidenced by the simulation of iKomissarov I 
i2004l) . Purely electrodynamic effects could not distort the 
field lines so that a flux tube that did not initially allow 
force-free inflow into the horizon could ever change. In this 
limit, in which the plasma never gains any substantial iner- 
tia, the field lines never get dragged inward into the hole by 
a transient. There was no evolution to force- free fiow into 
the black hole on these fiux tubes. 



5.3.2 GHM Field Configurations 

If the field lines do not allow an approximately force-free 
fiow into the horizon, the induced stresses imposed by the 
plasma creates a GHM dynamo in the ergopshere, by def- 
inition. Thus, the transients described above are transient 
GHM dynamos for example. The question posed by KDE is, 
are there ever any stable GHM dyi iamos? 

The simulation of fK omissarovj (|2004h and the analytic 
model of lPunslv and Coronitl 1.1990t) are examples of stable 
GHM field configurations. In these cases a vertical magnetic 
fiux through the equatorial plane prevents the plasma form 
fiowing towards the horizon. The field has to be stretched ra- 
dially, twisted azimuthally and possibly crossed by the ten- 
uous plasma in order for the fiow to proceed inwards. The 
struggle between the strong field and the relativistic iner- 
tia imparted to the plasma bj^the gravitational fiel d is the 
source of the GHM dynamo iPxmslv. (,200 ; Semenov et aJ 
J2004fl . 

5.3.3 Field Configuration: 2-D versus 3-D Time Evolution 

The scenario envisioned in ISorxiit and Uzdenskvl l|200 ,'j) re- 
quires 3-D. There needs to be an extra dimension to allow 
the degree of freedom that permits evolution of the mag- 
netic fiux tubes somewhat independent of the mass accre- 
tion. The flux tubes need interchange instabilities in order 
to allow plasma to accrete around them (i.e., so they can 
"swim" in the accretion flow). In 2-D ideal MHD, the ac- 
creting plasma will always build up in density until it can 
push the fleld lines ahead of it by ram pressure or just drag 
the field lines threading it into the horizon. This process will 
continue until the magnetic pressure roughly equals the in- 
going ram pressure. The plasma can never go around a flux 
tube. Conceivably, this extra degree of accretion freedom in 
3-D can produce a different fiux configuration for scenarios 
such as in Spruit and Uzdcnskv (2005); Dc Villicrs (2005). 

Perhaps this is what is happening to some extent in 
KDE. It is a much longer ergopsheric path to the horizon 
in terms of proper distance when a/M — 0.998 than the 
KDP simulation, a/M = 0.9 (a factor of 2.2) and the in-fall 
time to the asymptotic space-time near the horizon is there- 
fore much longer. Initially, the funnel might fill up as in the 
a — 0.9M case. However, late arriving flux tubes might have 
lost so much mass density to the outflow before they reach 
the asymptotic space-time near the horizon that they might 



actually attain a tenuous plasma state (i.e., Ua » 1) be- 
fore they arrive. As such, there is never enough inertial force 
to completely preconditioned the flux tube for future tenu- 
ous approximately force-free inflow. In GHM, the inertial 
dynamo forces are generated as the field torques the plasma 
on the way to the horizon since the azimuthal twisting is in- 
sufficient in the preconditioned field t o allow approxi mately 
force-free infiow lsemenov et aJ j2004ll : |Pun'sivl ^|2001^ . Even 
if there is an outfiow on these flux tubes, further accretion 
can still proceed in 3-D by moving around the flux tubes 
and a signiflcant GHM effect could coexist with the high 
accretion rate the deflnes quasar activity. 



5.4 Non Radial Field Geometry 

We verify the near radial fleld geometry in the funnel in 
figure 2 which is a plot of the poloidal field direction aver- 
aged over (f) in a, constant time snapshot. The poloidal field 
is radial on the boundary surfaces 2 and 4 of the Gaus- 
sian pillbox of figure 1. This does not support the notion of 
flowing along the fleld lines from the funnel wall, as dis- 
cussed in section 3.3, then following abrupt bends of the fleld 
lines toward the radial direction, so as to mimic a source of 
radial Poynting flux at 1.2M < r < 1.4M. Near the hori- 
zon, at lower latitudes in the funnel, 50° < 9 < 55°, where 
the strongest source is located (the red area in figure 
1), the field is almost exactly radial all the way to the in- 
ner boundary. A small deviation from a radial poloidal field 
near the horizon would not be able to produce strong source 
of in Poynting's theorem. Furthermore, if the poloidal 
field randomly oscillates about the radial direction (due to 
turbulence or numerical artifacts), near the horizon, then 
the oscillations in the Poynting fiux would average to zero 
in the time averaged plot in figure 1. At higher latitudes 
in the funnel, 30° < 9 < 50°, there is a region of rapid 
poloidal field variation right near the inner boundary. The 
fiuid is clearly highly turbulent, possibly a sign of reflections 
from the inner computational boundary. If the strong tur- 
bulence is physical, combined with the Ua >> 1 condition, 
this would provide strong magnetic pressure gradients that 
prevent the tenuous plasma from flowing force-free into the 
horizon. This is the condition described in section 5.3.2 that 
is most conducive to a GHM dynamo. 



6 CONCLUSION 

In this paper, we studied simulations of a magnetically dom- 
inated funnel of a rapidly rotating black hole, a/Af — 0.998. 
A source region that is responsible for creating over 72% of 
transported through funnel during the interval 2000A'/ < 
t < 8080 Af of the simulation was identified. Similarly, one 
can conclude that < 28% of emerging from the ergo- 
spheric funnel is from an inner boundary source, near the 
horizon. The small residual injected from the boundary 
into the accretion wind can be considered of electrodynamic 
origin. The distribution of in figure 1 is in contrast to the 
Blandford-Znajek solution in which essentially all the is 
of electrodynamic origin, i.e., it emanates from the horizon 
and passes through the accretion wind with minimal inter- 
action, thereby maintaining a virtually constant value along 
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Figure 2. A time snapshot of the averaged poloidal magnetic field direction. The arrows indicate direction not magnitude. The 
region is the same as that which is plotted in figure 1. The red line indicates the approximate boundary of the funnel. The arrows end 
approximately on a circle that is the inner calculational boundary. The plot is provided courtesy of John Hawley. 



each poloidal fiux tube throughout the ergosphere. The ex- 
istence of the electromagnetic source within the simulation 
was corroborated by existence of a cospatial electromagnetic 
angular momentum source using an entirely different data 
reduction method. Bo th of these results f ollow directly from 
the numerical data in lKrolik et all ((200^. 

There were two possible physical sources for in the 
funnel, there is the J • E source (associated with GHM in- 
ertial currents) and . By time averaging Faraday's law it 
was shown that |< F^e >!3> M |< F^t >|- The only source 
for both and < > in Poynting's theorem that is con- 
sistent with Faraday,s law was shown to be J • E = FtgJ^ . 
It was also pointed out that the decay of stored field en- 
ergy, de/dt could be the primary source of . This would 
suggest that the simulation is spontaneously creating energy 
and might be fatally flawed. In general, there could be a mix 
of numerical artifacts, and J-E that contribute to Poynt- 
ing's theorem. This can not be sorted out without further 
numerical work. 

The analysis presented above suggests some interest- 
ing physical possibilities and highlights the need for fur- 
ther analysis of 3-D simulations about rapidly rotating black 
holes (a/M = 0.998) in order to clarify the physics that cre- 
ates . Ideally, the 3-D simulation should be redone in 
ingoing Kerr-Schild coordinates (to eliminate boundary re- 
flections a the horizon) and preferably a conservative code. 
At present, there is no code efficient enough to handle the 
large computing times required. Whether this approach or 
theJDc Villicrs and Hawlcv (2003) approach is used, at least 
three consecutive time snapshots are needed in order to find 
the current distribution from Ampere's Law at each coarse 
time step data dump. It is likely that is highly turbu- 
lent near the event horizon, so it would be important to 



then time average the current distribution. Higher resolution 
(particularly in the 9 direction) might be required to resolve 
the issue of whether the source in figure 1 is from J^E±, 
or dt/dt. It would be prudent, if further calculations in 
Boyer-Lindquist coordinates are pursued, to locate the in- 
ner calculational boundary inside of the fast magneto-sonic 
surface at r « 1.065M (inner order to overcome potential 
reflections). Such future detailed analysis would determine if 
there is some new interesting physics or just some unforseen 
numerical errors in KDE. 



On a more general note, even though the simula tions of 
iDe Villiers et_a l (2003); Hirosc ct al (2004); Dc Villi ers et a 
|2005a"b'); Krolik ct al (120051) : iMcKinnev and Gammie 
|2004J ; .McKinnev . (.20051) are very impressive mathemati- 
cally for their ability to self consistently solve a very compli- 
cated problem, one must ask if they are representative of the 
astrophysical world. It was discussed in ^cKinncv ( 200^ 
that a radial field builds up in the funnel for any seed field 
configuration in the torus except for a purely toroidal field. 
This implies something generic about the result. The radial 
field always grows to seek a pressure balance with the infiow 
at the inner edge of the disk, therefore the mass accretion 
rate sets the field strength and the power output of the fun- 
nel. If this is a generic property of accretion disks then why 
are the vast majority (~ 90%) of quasars radio quiet. This 
tight coupling of the accretion rate to jet power seems to be 
a serious drawback of the results. From an astrophysical per- 
spective, it would seem prefera ble to develop simulat i ons o f 
scenarios like those described in lSoruit and Uzdenskvl (|200^ 
in order to describe the full panoply of quasar radio emis- 
sivity properties. 
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